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Abstract: Fundamental concepts in the quasi-one-dimensional geometry 
of disordered wires and random waveguides in which ideas of scaling and 
the transmission matrix were first introduced are reviewed. We discuss 
the use of the transmission matrix to describe the scaling, fluctuations, 
delay time, density of states, and control of waves propagating through 
and within disordered systems. Microwave measurements, random matrix 
theory calculations, and computer simulations are employed to study the 
statistics of transmission and focusing in single samples and the scaling of 
the probability distribution of transmission and transmittance in random 
ensembles. Finally, we explore the disposition of the energy density of 
transmission eigenchannels inside random media. 
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1. Introduction 

The scattering of waves by disorder is more the rule than the exception in condensed mate¬ 
rials. Scattering is the source of electrical resistance and makes materials opaque to classical 
waves. Whereas, all ordered samples of a given material are alike, every disordered sample is 
disordered in its own way. This can be seen directly in the distinct grainy pattern of intensity of 
scattered laser light for each sample 0. The speckle pattern within or beyond an illuminated 
sample arises from the interference of waves following all possible paths to any point. Though 
the wave field appears to be completely random in disordered samples, it is still possible to 
manipulate the incident wavefront to control the net transmission and to focus the transmitted 
wave. This was first achieved by focusing light to a point utilizing feedback from the trans¬ 
mitted wave to manipulate the input wave 0. In principle, a considerable measure of control 
over the output wave can be exerted by tailoring the incident wave if the field transmission 
coefficients between all incident and outgoing channels are known [3j. These coefficients are 
the elements of the transmission matrix (TM). In this review, we consider the use of the TM to 
describe the propagation and control of classical waves in random media. This will be placed in 
the electronic context in which the TM in the quasi-one-dimensional (quasi-ID) wire geometry 
was first utilized by Fisher and Lee Q to describe electronic conductance and by Dorokhov to 
describe the scaling of the eigenvalues of the transmission matrix and the conductance |5][6j. 
Because of the strong analogy between electronic and classical wave transport, it is possible to 
demonstrate fundamental predictions for electron transport in equivalent settings for classical 
waves in cases in which the predictions have not been demonstrated for electrons. In addition, 
consideration of the TM for classical waves has stimulated the investigation of universal aspects 
of transmission eigenvalues such as the dynamics and spatial profiles of the energy density of 
transmission eigenchannels, which cannot be explored directly in conductance measurements. 

The TM is a powerful tool for studying transport in disordered media in which the random 
field pattern excited by a monochromatic wave is coherent in time throughout the sample |^|9j. 
For classical waves, thermal fluctuations in position of the sample’s constituents in static sam¬ 
ples are much smaller than the wavelength so that the wave may be temporally coherent even 
on a macroscopic scale. In contrast, the electron wave is only phase-coherent in micron-sized 



conductors at ultralow temperature. Such samples lie between the microscopic atomic scale and 
macroscopic lengths. The suppression of diffusion m and the enhancement of fluctuations in 
conductance [T ZHH) due to interference in samples of intermediate or mesoscopic dimensions 
have been the focus of studies of electronic transport. 

Anderson 1101 showed that electrons in an unbounded lattice would be exponentially local¬ 
ized when the range of random values of the electron energy at a site is sufficiently large even 
when they are not trapped energetically. This leads to the absence of electron diffusion. Unlike 
the speckle pattern that is washed out by averaging over disorder, Anderson localization is an 
interference effect that survives averaging over random configurations. The suppression of av¬ 
erage transport is a consequence of the constructive interference of partial waves that follow the 
same path returning to a coherence volume in the sample but in time-reversed order 114 151, 
as illustrated in Fig. 1. The transverse dimensions of the coherence volume is the field correla¬ 
tion length A/2. Because the amplitudes and phases associated with such time-reversed partial 
waves are identical, the amplitude for return for paths following the loop for both senses is twice 
that for a single partial wave following the loop in either direction. The probability of return 
of a particle is therefore four times that for a single path and thus twice the probability given 
by the incoherent sum of probabilities for the two partial waves of each loop. This is the basis 
of weak localization [7j which reduces average transport. When the probability of return to a 
coherence volume reaches unity, the wave becomes localized. The effect of weak localization is 
observed directly as coherent backscattering in the domain of the reflected wavevector 116 -181. 
Reflection is enhanced by a factor of two over the background level in the retroreflection di¬ 
rection. The angular distribution of backscattered light is the Fourier transform of the intensity 
distribution of the point spread function for light incident on the incident surface. The enhance¬ 
ment falls in an angular range of 1 /k£, where k is the wavenumber and i the transport mean 
free path. Measurements of coherent backscattering therefore allow for the determination of the 
mean free path. 



Fig. 1. Illustration of a wave returning to a coherent volume of (A/2) 3 in three dimensions 
and its time-reversed partner. 


Ioffe and Regel [T9j suggested that, an electron cannot be considered to be freely propa¬ 
gating if the scattering mean free path were only a fraction of the wavelength. This yields the 
localization criterion in three dimensions, k£ < 1. This condition is also given by a calculation 
of the return probability based on diffusion theory in three dimensions. Subsequently, Thou- 
less [20 2T) argued that localization in any dimension should depend only upon the dimen¬ 
sionless ratio of the average width and spacing of energy levels, which has come to be known 
as the Thouless number, 8 = SE/AE. Electronic levels correspond to quasi-normal modes or 
resonances for classical waves in open systems. The width of levels, which in units of angular 














frequency is equal to the leakage rate of energy from the sample, is closely linked to the sen¬ 
sitivity of level energies to changes at the boundary since both are proportional to the ratio of 
the strengths of the mode at the boundary relative to the interior of the sample. When the wave 
is localized within the interior of a sample, the amplitude of the wave is exponentially small 
at the boundary and the mode is weakly coupled to the surrounding medium. Its lifetime is 
then long and its linewidth correspondingly narrow, so that, 8E < A E. On the other hand, when 
the wave is diffusive, modes extend throughout the medium; the wave then couples readily to 
its surroundings and the width of a typical level overlaps several modes, SE > A E. Thus, the 
localization threshold occurs at 8 = 1. 

Thouless 120 23 showed that for diffusive samples, 8 is equal to the conductance in units 
of the quantum of conductance 8 = g= G/(e 2 /h). Since the level width falls inversely with 
the square of sample length for diffusive samples while the spacing between levels falls only 
inversely with the sample length in the quasi-ID wire geometry of constant cross section, 8 falls 
inversely with sample length. Thus there will always be a crossover from diffusive to localized 
wave transport as the length of a wire increases. The crossover will occur at a sample length 
L = tq = AT, at which 8 = g = 1 [ 21], Here, tq is the localization length of the sample and N 
is the number of channels. The channels can be the propagating transverse modes in the ideal 
leads attached to a conductor at a given voltage. The total number of channels is approximately 
2kA /A 2 including two electronic spin states, where A is the constant cross sectional area of 
the sample and A is the electron wavelength. Abrahams et al. {22) argued that the variation 
of average g with the dimensions of the sample should depend only upon the value of g and 
the dimensionality. A localization transition is predicted to exist only in 3D in which g falls 
essentially exponentially with increasing sample size for g < 1 and increases for g> 1 {22) . 

Since the conductance exhibits large fluctuations for localized waves, the mean value of 
the conductance is insufficient to characterize the nature of electronic transport and the full 
distribution of the conductance needs to be given. Anderson et al. | 231 hypothesized that the 
logarithm of the conductance of localized waves in one-dimensional (ID) samples follows a 
Gaussian distribution with a variance approaching the average of the logarithm of the conduc¬ 
tance, (lng) = -L/q. The scaling of the distribution of conductance for localized waves in ID 
would then be completely specified by the single parameter L/B,. 

Electronic conductance in disordered ID conductors and transmission of classical waves 
in mesoscopic media are seen to have equivalent descriptions in the ID Landauer relation, 
g = T, where T is the transmission 124 [25) . The Landauer relation was extended to multi¬ 
channel mesoscopic systems via the TM {4)7 The elements f/ M of the transmission matrix t are 
the field transmission coefficients between incoming channels it a and output channels b. The 
dimensionless conductance in quasi-ID samples is equal to the sum of N 2 pairs of intensity 
transmission coefficient, which is the transmittance, T. g = b T«| 2 — T. The dimensionless 
conductance or transmittance can also be expressed as the sum of the transmission eigenvalues, 
g = Ti- Here, the T„ are the eigenvalues of the matrix product tt'. The transmission eigen¬ 
values can also be obtained via the singular value decomposition of the TM, t = UAV' , where 
U and V are unitary matrices with columns which are the singular vectors yielding the output 
and input of the transmission eigenchannels and A is the diagonal matrix with elements y/Zn. 

Dorokhov m considered the scaling of g in samples made up of N parallel ID disordered 
chains with transverse coupling to create a quasi-ID sample. He found that, instead of a sin¬ 
gle localization length !q, there exists N auxiliary localization lengths, !q n , which describe the 
scaling of the averages of each of the transmission eigenvalues vs. length via the relation, 

T„ = l/cosh 2 (L/4,)- (1) 


The inverse auxiliary localization lengths are are given by 1/<!;,, = (In — 1) /2q for integers 
n < A/2 so that these are equally spaced with spacing \/<q n+ \ — 1 jBq n = \/£, |5j|6|26-30|. 








An intuitive Coulomb gas model 1 26f3l) was developed to describe the statistics of transmis¬ 
sion eigenvalues. This is based on the model originally introduced by Dyson [32) to visualize 
the logarithmic repulsion between eigenvalue of a large random Hamiltonian matrix [33) . In 
the Coulomb gas model, the x n are associated with the positions of parallel lines of charges 
at x n = L/t^n = (2 n — l)L/2^) for n < N/ 2 and image charges of the same sign at —x n 1311 . 
The line charges in the model are embedded in a compensating continuous charge background 
of opposite sign. The background charge tends to screen the line charges from the interaction 
with other line charges. For diffusive waves, the spacing between neighboring charges is small 
compared to the screening length. As a result, there exists a logarithmic repulsion between the 
transmission eigenvalues due to the Coulomb interaction. The roughly uniform distribution of 
x„ imposed by the mutual repulsion of charges leads to the bimodal distribution of transmission 
eigenvalues x with peaks in the probability distribution for x near unity and at exponentially 


small values of T, p(x) = 
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. The repulsion between charges for diffu- 
Huctuitions. In quasi-ID disordered systems with 


sive waves leads to universal conductance 
time-reversal symmetry, var(T) =2/15 [28]. 

Universal conductance fluctuations were first measured in mesoscopic conductors at low 
temperatures HMD- Fluctuations of conductance of a gold wire at 10 mK as a function of 
magnetic field applied perpendicular to the wire are shown in Fig. 2. Dorokhov |6J showed 
that transport through disordered conductors is dominated by g eigenvalues with > L and 
so transmission through these channels is of order unity. The number of such open eigenchan- 
nels with x n > 1/e, which is essentially the conductance for diffusive quasi-ID samples, is 
inversely proportional to the sample length in accord with Ohm’s law. For localized waves, 
however, the spacing between charges in the Coulomb gas model is greater than the screening 
length so that the Coulomb interaction between charges is diminished. The smallest of the x„, 
x\, is larger than unity for localized waves so the transmission eigenvalues fall exponentially 
with n, (x n +i)/(x n ) = exp(— 2L/t j ). As a result, the conductance is dominated by the highest 
transmission eigenchannel X\ [6|[26j[27 29 -31). 
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Fig. 2. Fluctuations of conductance of a gold wire at 10 mK as a function of magnetic field 
applied perpendicular to the wire. The wire is of length 310 nm and width 25 nm. (From 
Ref. (37)) 

In contrast to electronic systems, for which only the conductance can be measured, it is in 
principle possible to measure the elements of the TM for classical waves. Thus the transmission 
coefficients for the field, tb a , intensity Tf >a = \tba\ 2 , total transmission T a = Tba and transmit¬ 
tance T can be measured. Spectra of these variables for a single sample supporting localized 
and diffusive waves in different frequency ranges in which g =6.9 and g = 0.37 are shown in 



























Fig. 3. The statistics of random ensembles of statistically equivalent samples can be obtained in 
microwave measurements for a succession of sample realizations produced by briefly rotating 
a copper tube containing randomly positioned dielectric spheres. Ensembles of random sys¬ 
tem can be studied optically by illuminating different regions of a random slab. The statistics 
of internal motion of colloidal samples can also be explored in measurements of the temporal 
correlation of transmitted light [38} 40). 
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Fig. 3. Spectra of the transmission coefficients normalized to ensemble average values. 
In- and out-of-phase field transmission coefficients, field E^/\fJTboJ), transmission co¬ 
efficients of intensity Sb a = Tj, a /(Ti, a ), total transmission s a = T a /{T a ) and transmittance 
s = T/ (T) for microwave radiation propagating through a random waveguides for a diffu¬ 
sive sample (left column) and a localized sample (right column). For diffusive waves, 8 > 1, 
many modes contribute to transmission at all frequencies and for all source and detector 
positions. In contrast, for localized waves 8 < 1, distinct peaks appear when the incident 
radiation is on resonance with a quasi-normal mode. The resonance condition holds for all 
source and detector positions and therefore sharp peaks remain even when transmission is 
integrated over space. 


The first question that was addressed in bounded samples was the scaling of conductance. 
This is practically difficult to realize in electronic measurements because of the limited scale 
of mesoscopic samples and the difficulty of producing statistically equivalent samples. The 
scaling of optical transmission in the slab geometry can be measured in the slab geometry by 
forming a wedge of random material and focusing a laser beam on the wedge as it is translated 
perpendicular to its vertex ]4l]|42) . The transmission of light through a slab of rutile titania 
particles in a polystyrene matrix is seen in Fig. 4. The lower envelopes of the transmission 
curves in Figs. 4a,b fall inversely with sample thickness L up to a length equal to the absorption 
length L a indicating that the wave is diffusive. The absorption length is given by L a = \/Dx a , 
where the diffusion coefficient is D = v e £/3, v e is the energy velocity in the medium and z a is 
the absorption time. Spikes in transmission occur at spots in the sample at which a pit has been 
gouged out of the otherwise smooth surface in the polishing process. Spectra of transmitted 

















intensity in the far field at different sample thicknesses are shown in Fig. 4(c). The correlation 
frequencies of fluctuations in transmission at different sample thicknesses give the average 
linewidth of electromagnetic modes of the medium 1431. In this diffusive sample, the correlation 
frequency scales as L 2 for L<L a |43] |. The field correlation function and the temporal profile 
of transmission form a Fourier transform pair and give the diffusion coefficient and absorption 
length which is consistent with measurement of the scaling of static transmission 63 - 



Fig. 4. Scaling of optical transmission through a wedged sample consisting of TiCL em¬ 
bedded in polystyrene. (a,b) Log-log and semi-log plots of the scaling of the transmission 
coefficient, T(L). The inverse variation of T(L) with L is seen in the log-log plot and the 
exponentially falloff of transmission beyond the absorption length L a is seen in the semi¬ 
log plot. The absorption length L a = 112 ±5 fim of the sample is indicated by the dashed 
line in both figures. The solid line is a fit of the expression T(L) = 5a£)/vsinh(aL) to 
the data. Here D is the diffusion coefficient, v is the speed of light in the medium and a 
is the inverse of L a . (c) Spectra of transmission for various sample lengths. The measured 
transmission in these experiments never reaches zero because the polarized intensity only 
vanishes at point singularities while measurements are made with a finite aperture placed 
before the photomultiplier tube. (From Ref. ED) 


Establishing that a wave is localized through scaling is a challenge in the presence of absorp- 

481. But the 


tion of classical waves and of dephasing of the electronic wave function 142 441 


magnitude of fluctuations of intensity and total transmission normalized by their ensemble aver¬ 
age values give a straightforward measure of the degree of localization which is robust against 
absorption |48) . The magnitudes of relative fluctuations in electronic conductance 6363 - 
classical intensity and total transmission and the degrees of spatial, spectral and temporal cor¬ 
relation of flux increase through the localization transition 140 43 49 63], Large fluctuations 
in transmitted intensity were observed in the critical regime of the Anderson localization tran¬ 
sition in acoustic measurements in a slab of brazed aluminum beads [:64j. These results are 
consistent with the agreement found between the non-exponential decay of a transmitted pulse 
measured in this sample and self-consistent localization theory |65j|. Microwave measurements 
in quasi-ID samples of the time-dependence of transmission in more deeply localized sam¬ 
ples |66j show dramatically reduced decay rates for transmission compared to both diffusion 
theory and the self-consistent localization theory. Since the slow-down of decay rates is as¬ 
sociated with long-lived modes |67 69), this indicates that the self-consistent theory is valid 
near the localization transition but not for deeply localized waves 165 70-72). The approach 
to localization can also be seen in the saturation of the spread of a wave with time delay on 
the output surface of the sample, which is independent of absorption. This was observed in 
acoustic measurements |64| and in pulsed laser experiments in random slabs of rutile titania 
particles 63 - The local strength of scattering in the optical experiments was determined from 
measurements of the width of the coherent backscattering cone 173174) . Values of the product of 



























the wave vector and the transport mean free path, kt were larger than the value predicted by the 
Ioffe-Regel criterion 1191 for localization in three dimensions of M < 1. The lowest value of k£ 
for these samples was 2.5. The Ioffe-Regel criterion is consistent with the Thouless criterion for 
localization |2j 751. Localization in samples with kt> 1 might occur in random samples with 
residual order |76j or in samples with resonances in individual scattering elements 77 -80). In 


such samples, the density of states (DOS) may fall below the uniform medium value 8 nv 2 /Vp, 
where v is the frequency, v p is the phase velocity. Localization might then occur from values 
of k£>\. 

Additional evidence that these samples are in the critical regime is seen in the continuous 
slowing decay of pulsed transmission in these samples. This shows that the exponential decay 
of transmission with time delay expected for diffusive samples breaks down at long times. 
Departures from an exponential decay of transmission for diffusive waves with values of g~ 2.5 
have been observed in microwave transmission (68). This indicates that the distribution of mode 
lifetimes increases even in diffusive samples as the crossover to localization is approached 
| 67]|68) . An important question is whether the time dependence of the transverse spatial profile 
represent the reduced spatial extent of some pre-localized modes as the localization transition is 
approached or localization 167 81 82]. Transverse optical localization |83] has been observed 
in samples that are uniform in the longitudinal direction but disordered in two [:84] or one 
transverse dimension. Whether it is possible to localize electromagnetic radiation in three 


dimensions in samples without residual order remains an open question [76.86). 

The TM has been recently measured for light 1 3|87 89), microwave radiation J29) and acous¬ 
tic waves [90)91] in static samples. The TM gives a full description of the degree to which the 
incident wavefront can be manipulated to control the transmitted wave ]3)[92). This may be ex¬ 
ploited to focus radiation and enhance or suppress transmission for applications in imaging and 
telecommunications 1 3|29|87 - 89)93 - TT2] . The ability to control propagation of classical waves 
in a random medium was first demonstrated in acoustics by focusing a pulse within and through 
a scattering medium 013 • This is achieved by time-reversing the signal recorded by an array 
of transducers. The measured time-varying signal is played back in time and a pulse emerges at 
the location of the source. For monochromatic radiation, it is possible to focus the transmitted 
wave by shaping the incident waveform. Vellekoop and Mosk @ demonstrated optical focus¬ 
ing through an opaque sample by adjusting the incident wavefront reflected from a spatial light 
modulator using a genetic algorithm with feedback from the intensity at the focus. When the 
TM of the scattering sample is determined, it is possible to focus at any desired channel on the 
output surface. This can be achieved by producing an incident waveform tp a , which is the phase 
conjugate of the field one would obtain on the incident plane when a source is placed at a point 
/I on the output plane. This would bring the fields from different incident channels a arriving 
at ji in phase so that the intensity at the point is greatly increased by constructive interference. 
Vellekoop and Mosk [ 921 found that the background intensity in the transmitted wavefront for 
optimally focused radiation is also enhanced relative to the average transmitted intensity. In ad¬ 
dition, the total transmission can be significantly enhanced to give nearly complete transmission 
when the incident wavefront corresponds to the highest transmission eigenchannel. 

The TM gives the fullest account of transmission and conductance and enables the control of 
the transmitted flux. However, the TM cannot provide the intensity profiles of the transmission 
eigenchannels inside the sample. Recent simulations by Choi et al. |1 14) explored the spatial 
profile of the intensity of eigenchannels inside a single random sample. They have observed 
that the energy stored within high transmission eigenchannels is enhanced. This was confirmed 
by Gerardin et. al. (ED in measurements of flexural waves propagating in a disordered elastic 
waveguide. Using a laser source and a heterodyne interferometer, they have measured the full 
scattering matrix of the medium. After a proper normalization of the measured scattering ma- 































trix, they found that the transmission eigenvalues follow a bimodal distribution in agreement 
with theoretical prediction. They have also demonstrated that the intensity associated with high 
transmitting eigenchannels increases within the medium with a peak located near the middle of 
the sample. This was presented in Fig. 5. We have shown recently in microwave measurements 
that the integral of the intensity inside the sample for a specific eigenchannel yields the dwell 
time of the eigenchannel [115). The dwell time of the each of the eigenchannel is also equal to 
its contribution to the DOS of the sample so that the sum of the dwell times for N eigenchannels 
gives the DOS of the disordered material. 



Fig. 5. Absolute value of the wave field in the scattering medium at f = 0.36 MHz associated 
with (a), an incident plane wave, (b), a closed eigenchannel, (c), an open eigenchannel 
deduced from the measured S matrix, and (d), an open eigenchannel deduced from the 
normalized matrix S. The corresponding intensities averaged over the wave guide section 
(y axis) are shown versus depth x in lower panels (e)(h). They are all normalized by the 
intensity at the plane of sources (x = 0). (From Ref. HD) 


The average intensity falls linearly inside the medium in diffusive samples, as required by 
the Fick’s law. For a localized sample, wave interference cannot be neglected in calculating the 
average profile ]116|117| . The intensity profile for localized waves excited by radiation incident 
from one side of the sample is found to fall more rapidly toward the center of the open medium 
170 - 72) . This profile can be found by solving a generalized diffusion equation with a position 
dependent diffusion coefficient, which reflects the increasing renormalization of transport with 
depth into the sample due to wave interference Cg. Recent simulations of intensity profiles 
inside the material for a random ensemble have suggested a universal expression for the profile 
of the energy density of the transmission eigenchannels 11181. These profiles may be expressed 
in terms of the auxiliary localization lengths of the corresponding transmission eigenchannels. 
The structure of these profiles is given by solutions of the generalized diffusion equation [1181. 

In the sections of this review that follow, we will focus on recent microwave measurements 
and computer simulations of spectra of the TM in quasi-ID random waveguides in the crossover 
from diffusive to localized waves. In Section II, we discuss the statistics of transmission in sin¬ 
gle sample realizations and the statistics of the transmittance in ensembles of random samples. 
The statistics in singe samples with given transmittance T are found to depend upon only a sin¬ 
gle parameter, which is the eigenchannel participation number, M = (Yn t «) 2 /L^=i t ,7 1191. 






























The probability distribution of T is observed to change from a one-sided log-normal distribu¬ 
tion to a log-normal distribution as M -> 1 as the wave becomes deeply localized. In addition, 
the statistics of transmission are found to approach predictions for ID. In particular, the single 
parameter scaling (SPS) hypothesis predicted for ID random samples is found experimentally 
to hold for deeply localized waves in a quasi-ID geometry. Focusing of both static and pulsed 
radiation through and within opaque samples is discussed in Section III. Expressions for the 
spatial profile of the focused beam and the contrast between the focused and background in¬ 
tensities are demonstrated. In Section IV, we show that a complementary set of parameters to 
the transmission eigenvalues provides the dwell time and the contribution of each eigenchannel 
to the DOS. These parameters correspond to the integral of the energy density profile inside 
the sample for the corresponding transmission eigenchannel. A summary of recent results and 
prospects for applications are given in in section V. 

2. Statistics of transmission in single samples and ensembles of random samples 

The statistics of transmission in a single disordered sample are of particular interest since appli¬ 
cations are typically in a specific sample. Such statistics are an essential part of the statistics of 
transmission over a collection of random samples. The transmitted field with a given polariza¬ 
tion for an incident channel a in a single random speckle pattern is a complex Gaussian random 
variable. Since the transmission is the square of the field on the output surface, T bcl = \t ba \ 2 the 
corresponding probability distribution of intensity relative to the average in the speckle pattern, 
Tba/(Zb T ba/N) = NTba/Ta, is a negative exponential function, P(NT ba /T a ) = exp {-NT ba /T a ) 

1 120[|l2l| . Thus, the statistics of the transmitted flux relative to its average in a single TM are 
fully specified by the statistics of total transmission T a normalized by its mean within the TM, 
T a /(T a ) a = NT a /T. Here, (.. ) a indicates the average over all incident channels a in a single 
sample. 



Fig. 6. The variance of the normalized total transmission vs. the inverse of the eigenchannel 
participation number M is obtained from measurements by grouping single samples with 
the same value of A/ -1 from ensembles with a wide range of values of g. The solid line is 
var (NT a /T) =M~ l . (From Ref. (Tl9|) 

Using the singular value decomposition of the TM, the variance of the total transmission for 
an incident channel a relative to the average over all possible incident channels can be expressed 
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For diffusive samples without absorption, the bimodal distribution of the eigenvalues gives 
( M) = 3g/2. The value of g for diffusive waves is close to the number of open eigenchan- 
nels with transmission eigenvalues greater than 1/e [6 1221. For deeply localized waves, the 
number of open eigenchannels falls to zero, but M approaches unity since the transmittance is 
dominated by a single transmission eigenvalue. 

These calculations apply to a single large TM and to subsets of samples with the same values 
of M 1 . We consider measurements of samples with small value of N and group together sub¬ 
sets of samples with similar values of M . Figure 6 demonstrates that vai(NT a /T) is close to 
M _1 in subsets of samples with given M . 

The measurements of microwave transmission described in this review were carried out in 
collections of samples consisting of randomly positioned alumina spheres with refractive in¬ 
dex 3.14 within a copper tube. The dielectric spheres were embedded in Styrofoam shells to 
produce an alumina volume fraction of 0.068. Microwave radiation was launched and detected 
by wire antennas connected to a vector network analyzer. The orientation of the wire antennas 
determines the polarization of the field that is detected. The source and detector antennas are 
mounted on 2D translation stages and moved to positions on a square grid over the open ends of 
the waveguide at which spectra of the in- and out-of-phase components of the transmitted elec¬ 
tric field t bn are collected. Measurements are made over two frequency ranges of 10-10.24 GHz 
and 14.7-14.94 GHz. The lower frequency range is just above the first Mie resonance of the in¬ 
dividual alumina spheres and so the scattering is strong and the wave is localized in a relatively 
short length of the random waveguide. The wave is diffusive in the higher frequency range. The 
number of allowed propagating waveguide modes increases with frequency; they are ~30 and 
~66 in the centers of the two frequency ranges. Measurements of spectra of the TM in the two 
frequency ranges are made for sample lengths between 23 and 102 cm so that ensembles over a 
wide range of g are studied. The impact of absorption on the statistics of transmission is largely 
removed by Fourier transforming the field spectrum into the time domain and multiplying the 
time signal by exp('f/2T„), where t is the time delay and T a is the absorption time 146 48) . 

Measured distributions of relative total transmission and intensity for subsets of samples with 
M in the ranges 0.17±0.01 and 0.995 ± 0.005 are shown in Fig. 7. The distributions in sam¬ 
ples with the same values of M 1 but drawn from ensembles of different g are seen to coincide. 
This shows that the distributions P(NT a /T ) depend only upon the value of M. P(NT a /T) are 
also presented in Fig. 7(b) for M 1 near unity in measurements for two ensembles of localized 
samples with g = 0.37 and 0.17. When Ti T 2 , transmission is dominated by the first transmis¬ 
sion eigenchannel so that, t ba = £„ u bny /%,vj m ~ u bl s/T\v\ a and T ba = \u bl | 2 Ti |v{j 2 . This gives 
NT a /T = |via| 2 . The probability distribution P(NT a /T) is therefore expected to be a negative 
exponential, P(NT a /T) = e\p(—NT a /T). Measurements are in excellent agreement with this 
conjecture, as shown in Fig. 7(b). 

The eigenchannel participation number M can be linked directly to the degree of inten¬ 
sity correlation within a single sample. In the limit IV 1, infinite-range correlation due 
to fluctuations in T vanishes if computed for the TM of a single sample. The cumulant 
correlation function of transmitted intensity relative to its average, T //V 2 , is then CV , = 
([ TbaTb'a' ~ ( T/N 2 ) 2 ]/(T / N 2 ) 2 )m ■ The values of residual short-range and long-range correla¬ 
tion determined from the variances of relative intensity and total transmission are unity and 
M 1 , respectively, giving 


Cba,b'a' ^na’ &bb' ± Tf (&aa ! ± ) ■ 


(3) 









Fig. 7. Probability distributions of normalized total transmission and intensity in sub¬ 
ensembles of random samples with similar values of M _1 . (a) The probability distributions 
of normalized total transmission P(NT a /T ) for samples with M _l = 0.17 ±0.01 from two 
random ensembles with g = 3.9 (Green open circles) and g = 0.17 (red filled circles) are 
seen to overlap, confirming that P(NT a / T) depends only on the value of M. The solid curve 
is a calculation of P(T a /{T a )) with var(T a /(T a )) replaced with M _l in the expression from 
Ref. |59 l|60) . (b) P(NT a /T) for M 1 in the range of for two ensembles with g = 0.37 and 
0.17. The solid line is an exponential distribution, exp(—NT a /T). (c,d) The normalized 
intensity distributions, P(N 2 T/, a /T) for the same as in (a,b). The solid lines are the calcu¬ 
lations of the distribution based upon P(NT a /T) in (a,b) and the exponential falloff of the 
distribution of the normalized intensity P(NTi, a /T a ) = exp(—NTj, a /T a ). (From Ref. ) 119)) 


This leads to var(A ,2 7/ M /7) = 1 - 2M 1 . This is confirmed in measurements by the correspon¬ 
dence of the measured variances of relative intensity for the two values of M 1 of 0.17 and 
0.995 of 1.38 and 3.04 with the calculated values of 1.34 and 3.00. 

Since the statistics of transmission in a single sample at a given frequency depend only on 
M, they do not reflect the nature of transport within a disordered material. This would require 
knowledge of the distribution of M. Previous studies have shown that the crossover to Anderson 
localization can be tracked in terms of the fluctuation of the intensity and total transmission over 
collections of samples in an ensemble 148 58-60 62 63). In particular, for diffusive waves in 
the quasi-ID geometry, the variances of and T a normalized by their ensemble averages are 
linked to the dimensionless conductance g via the relation, var(sfo„ = 7* a /(7i fl )) = 1 +4/3 g 
and var(s« = T a /(T a )) =2/3 g, respectively. var(s a ) is also equal to the degree of long range 
intensity correlation, fc = where 8st, a = 1. In many circumstances, it is not 

possible to determine the absolute value of the ensemble average of the conductance, g, and 
in measurements in some geometries, such as optical measurements of transmission in a slab, 
the precise correspondence of transmission and conductance is not clear. Nonetheless, it is 
still possible to determine the degree of localization in a sample from measurements of the 
variances of normalized intensity or total transmission. Because of infinite-range mesoscopic 









































intensity correlation, the fluctuations do not self-average and remain prominent even for the 
most spatially averaged transmission quantity, the transmittance T fTTj{13||23)|40||123}{127| . 

Measurements of the TM for classical waves make it possible to explore long-standing pre¬ 
dictions regarding the statistics of conductance in mesoscopic systems in the crossover to lo¬ 
calization from the perspective of the transmission eigenvalues |23||128 -135]. Unlike measure¬ 
ments of the electronic conductance, for which the measurement of current gives the flux of 
electrons directly, the transmittance of classical waves is measured on a grid in real space. In 
such measurements, it is difficult to measure transmission for all incident and outgoing chan¬ 
nels. Goetschy and Stone Jl 36) have considered the impact of the incomplete measurements 
of the TM in momentum space upon the statistics of the transmission eigenvalues. They find 
that the distribution of the density of transmission eigenvalues changes progressively from the 
bimodal distribution for diffusive waves to a distribution characteristic of Gaussian random 
matrices |137| as the ratio of measured channels N' to the total number of channels N on the 
input and output sides, m\ = N[ /N and m 2 = N' 2 /N, decreases. We find the impact of incom¬ 
plete measurements upon the statistics of eigenvalues is low as long as N' > g [30 136 1381. 
Thus, the TM for localized waves determined in a single polarization in real space can faithfully 
represent the statistics of the transmission eigenvalues and the transmittance. 
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Fig. 8. Probability distribution of microwave conductance, (a) /’(lnT) for g = 0.37 (red 
dots) and 0.045 (green asterisk). The solid black line is a Gaussian fit to the data. For g = 

0.045, all of the data points are included, whereas for g = 0.37, only data to the left of the 
peak is used in the fit. (b) P(T) for g = 0.37 in a semi-log plot exhibits an exponential tail. 

(From Ref. |30}) 

In Fig. 8, we present the probability distributions of InT for microwave radiation transmitted 

























through random ensembles of dielectric samples with g = 0.37 and 0.045. Since the measured 
distributions of the normalized transmittance s = T/(T) depends only upon g for localized 
samples, the value of g for the localized samples can be obtained by comparing the probability 
distribution P(s) with numerical simulations for a given value of g. P(lnT) is seen in Fig. 8(a) 
to be close to a Gaussian distribution for g = 0.045. For g = 0.37, the low-transmission side 
of P(lnF) is seen to be well fit by a Gaussian distribution, while the high transmission side 
falls sharply above the peak at lnT = —0.5. Above T = 1.1, P(T) is seen in Fig. 8(b) to fall 
exponentially, in accord with predictions in |134) . We show below that in addition to explaining 
the origin of the distributions of transmittance in the diffusive and localized limits, the Coulomb 
gas model can explain the anomalous distribution of transmission in the crossover between 
diffusive and localized transport [26jj3l]|. Understanding transport in terms of the charge model 
is of further interest because the charge model provides a parallel treatment of diffusive and 
localized waves and so forms the basis for a universal description of wave transport. 



Fig. 9. Coulomb gas model of transmission eigenvalues and conductance, (a) Average posi¬ 
tions of charges and their images as the position of the first charge x\ change in the random 
ensemble with g = 0.37. The dashed lines show the average positions of the charges for 
this ensemble. The curly brackets in this case represents the average over a subset of trans¬ 
mission matrices with the specified value of x |. (b) Average positions of charges vs. lnT 
in the same ensemble. Here, the curly brackets indicate the averaging is over a subset of 
transmission matrices with the specified value of In T. (From Ref. |30| ) 

The variation of the average positions of the charges x n associated with t„ for different po- 
























sitions of the first charge xi in the random ensemble with g= 0.37 is plotted in Fig. 9. The 
repulsion between xi, which is associated with the highest transmission eigenvalue Ti, and its 
image at -x\ enforces a ceiling for T\ of unity. The average spacing between x\ and xi increases 
as the value of x\ decreases since the charge at xt then interacts strongly with the charge at x\ 
as well as with its nearby image at — x\. At the same time, the spacing between X 2 and X 3 and 
their average positions hardly change. This reflects the tendency to heal large fluctuations in 
charge positions for more remote charges. The source of the sharp cutoff in P(lnT) can be seen 
in the context of the Coulomb gas model by examining the spacings of charges for different 
values of In T, as shown in Fig. 9(b). A relatively high value of T is only achieved when the 
first charge is near the origin. This is an unlikely circumstance because this charge is strongly 
repelled by its image. P(\nT) would be expected to fall off especially rapidly for values of T 
above unity since this would requires two charges along with their images to be close together 
near the origin as seen in Fig. 9(b). 

The probability distribution of In T for the deeply localized waves is seen to follow a Gaus¬ 
sian distribution. We further explore the relation between the first two moments of the dis¬ 
tribution. The ratio of cr 2 = varfln T) and ■ (In T), vs. Ljq is plotted in Fig. 10 and seen to 
approach unity forL q ~ 24 cm. This occurs just as the weighted average of M 1 , (MT) / (T) 
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Fig. 10. Approach to single-parameter scaling in multi-channel random systems. Si = 
—var(ln7 1 )/(ln7’), {MT)/{T) and {M~ l T)/(T) are plotted vs. L/iq. The dashed line is 
the prediction of SPS in the limit of deeply localized samples in ID disordered materials. 
(From Ref. ffl) 


and (M _1 T) / (T) approach unity. The latter average closely tracks SSi. Thus, SPS is approached 
in a multichannel quasi-ID random system when transmission is through a single eigenchannel 
and transport is essentially one dimensional [23|[30||139||140| . In this limit, the value of M is 
close to one in nearly every realization of the random sample except for sample configurations 
with especially small values of T, and the relationship between the statistics of transmittance, 
total transmission and intensity is particularly straightforward. We can therefore obtain the 
distributions of total transmission and intensity from the distribution of the transmittance by 
setting M = 1 and P(T,M) = P(T). Because the distribution of normalized intensity for a given 
value of T a is a negative exponential, fluctuations of normalized intensity are closely linked to 
fluctuations of normalized total transmission, ( s n ba ) = w! (,v"). Similarly, in the limit of M = I, 
(s n a ) = n\(s n ), where s = T /( T). This gives the expression for the variance of st, a when trans¬ 
port is dominated by a single transmission eigenchannel, var(s/, a ) = 2var(s a ) + 1 = 4var(s) +3 










11411. We find the values of the variances of s ba , s a and s to be 27.5, 13.2 and 6.4 in the sample 
of L = 102 cm, which is consistent with these relations. This further confirms that the transport 
in this random ensemble is via a single eigenchannel. 

3. Controlling wave propagation in a scattering medium 

Since a random speckle pattern develops on the output plane of a multiply scattering sample 
illuminated with monochromatic light, it appears that the information regarding incident wave 
is completely washed out jT|. Nonetheless, because the transmitted field is linked to the incident 
field via the TM of the sample E b = tbaE a , it is possible to solve the inverse problem and to 
retrieve the incident waveform by recording the TM and the transmitted field j'4||6| |28]|34) . In 
particular, measuring the TM provides the specific waveform to create a focal spot |3j 88]|98 

ITOOliToH . 

Transmitted waves can be focused at a target point ji by applying an incident waveform 
which is the phase conjugation of the field transmission coefficients between the target point 
13 and the incident points a, t^ a 1921. A peak in the intensity emerges at /j, since the transmit¬ 
ted fields from different incident points a arrive at ji in phase and interfere constructively while 
these fields are randomly phased at other points. The incident field is normalized by |f/ 3 fl | 2 , 

so that the incident power is unity. The intensity at focal point ji is equal to YjJj] \tp a \ 2 < which 
is the total transmission Tp from the point ji to the opposite surface. This yields a factor of N 
enhancement over the ensemble average transmission (7j, a ). The background intensity is also 
increased because of the long range intensity correlation within the TM so that the contrast in 
focusing is smaller than N. The intensity at the focus and in the background can be written as, 


=I«T,|M„ /3 |-and4 /i3 = 


L, t n u nb u n p I 


, respectively 119 


ment of the TM is defined as, fj, = {Ip)/ (4), in which { 
background channels b for a given focusing point ji. In the limit of N 1, we find, 


. The contrast for a single measure- 
} b n indicates the average over the 


A* = 


1/M- 1/A' 


(4) 



Fig. 11. Contrast in optimal focusing vs. eigenchannel participation number M. The open 
circles and squares represent measurements from transmission matrices N = 30 and 66 
channels, respectively. The filled triangles give results for N' x N' matrices with N' = 30 
for points selected from a larger matrix of size N = 66. Phase conjugation is applied within 
the reduced matrix to achieve maximal focusing. Eq. 4 is represented by the solid red and 
dashed blue curves for N — 30 and 66, respectively. (From Ref. (HD) 












The results of measurements of the contrast in optimal focusing for diffusive samples of 
length L = 23 cm for incident waves in two orthogonal polarization are shown in Fig. 11 to be 
in excellent agreement with Eq. 4. When the number of measured points N' is smaller than N 
and therefore the corresponding eigenvalue participation number M' is smaller than M for the 
complete matrix, the contrast is given by Eq. 4 with the substitutions M —¥ M' and N —> N'. 
This is demonstrated in Fig. 11 with the contrast in samples with N = 66 but with the contrast 
computed only for N' = 30 points falling on the curve for N = 30. These results may be applied 
to optical measurements of the transmission matrix, where the full matrix is not determined. 

In contrast to a conventional focusing lens, in which the size of the focused beam is diffrac¬ 
tion limited and so equal to A /2 divided by the numerical aperture of the lens, the profile of 
the focused beam through a disordered medium depends only upon the property of the random 
system itself |56| 94 1011. The average intensity at a point b at a distance A r from the focal 
point (If oc (Ar)} normalized by the peak intensity for a diffusive sample can be expressed in 
terms of degree of the long range intensity correlation K and the square of the amplitude of the 
field correlation function, F(Ar) = \E(r)E*(r+Ar)\/(I(r))(I(r + Ar)}, 


(If oc {Ar))p _ F(Ar) + K 

(Tp) 1 + K ' 

The ratio between the background intensity F(Ar A) and the focus intensity F(Ar = 0) is 
equal to k/{\ + K"), since F(Ar) —> 0 for Ar 3> A. In the diffusive limit, K ~ 1 /M <C 1, this gives 
a contrast of fx = M and the focused profile is reduced to F(Ar). Since the intensity correlation 
function in the diffusive limit is F(Ar) according to the field factorization approximation, the 
profile is then expected to be the same as the intensity correlation, which has been demonstrated 
in optical measurements (94). 
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Fig. 12. The ensemble average of normalized intensity for focused radiation (blue circles) 
is compared to Eq. 5 (blue solid line) for L = 61 cm for diffusive and localized waves. 
The transmitted field for a single polarization is measured along a line with a spacing 
of 2 mm and the source antenna is translated along 49 points on a square grid with 9 
mm spacing. F(Ar) (blue dots) is fit with the theoretical expression obtained from the 
Fourier transform of the specific intensity (dashed blue line). The black dashed line is 
proportional to (7), a )/ (If oc (Ar = 0)) = l/N. Equation 5 is not valid for localized waves, 
but good agreement is obtained when K is replaced by 1 / (ju — 1) in Eq. 5 with /r determined 
experimentally. (From Ref. 


The normalized profile of the focusing beam via phase conjugation averaged over all possible 



















focusing spot and over all frequency points is presented in Fig. 12. The profile of the focused 
beam is in excellent agreement with Eq. 5 for diffusive samples. The black dashed line is 1 /N, 
which is the ratio between the ensemble averaged intensity and the intensity at the focus. The 
background intensity in optimal focusing is seen to be much greater than the average transmis¬ 
sion over a random ensemble and is enhanced by a factor of N/M. For localized waves, only a 
single transmission eigenchannel contributes so that the transmitted wave cannot be focused by 
applying phase conjugation. 

It is also possible to focus a wave inside a random medium. Since the field inside a material 
is in general not accessible, many virtual guide stars have been created inside the sample in 
order to facilitate the focusing at a desired location. These include the use of feedback from 
florescent particles embedded in the sample 195 |, photoacoustic |106| and acousto-optical ef¬ 
fects J97]|104| and the use of many nonlinear phenomena, such second harmonic generation. 
The contrast of focusing inside a disordered sample has been explored in simulations of a scalar 
wave propagating through Q1D samples which are locally two dimensional. The wave equation 
V 2 £’(x,y) + k^e{x,y)E{x,y) = 0 is discretized on a square grid and solved using the recursive 
Greens function method 11421. Here, ko is the wave number in the vacuum and e(x,y) is the 
spatial varying dielectric function. e(x,y) is set to be unity in the two ideal leads and is equal to 
1 ± 8s in the disordered sample. 8s is drawn from a rectangle distribution whose width repre¬ 
sents the strength of the disorder. Both internal and external reflection are negligible since the 
sample is index matched to the leads. The product of the wave number in the leads ko and the 
grid spacing is unity and the dimension of the samples are measured in units of the grid spacing. 
The Green’s function G(r, r) is calculated between grid points r = (0, v) on the incident plane 
and r' = (x , ._v / ) at depth x' within the interior of the sample with y the transverse coordinate. 
We considered an ensemble of diffusive samples with length L = 800 and width L, = 500. The 
fluctuation of dielectric function at each site is drawn randomly from a rectangular distribution 
from 0.8 to 1.2. This gives g = 16. 

Optimal focusing at a location at (xo,yo) inside the sample can be obtained by phase conju¬ 
gating the Green’s function between arrays of points in the incident plane and the target spot 
(xn. yo). In Fig. 13(a), we present the intensity profile of the focused wave at the center of the 
transverse dimension L t /2 of the waveguide at two depths x = L/4 and x = L/2. In both cases, 
intensity peak emerges at the focal spot with a contrast of ~ 37 and ~ 17, respectively. 




Fig. 13. Optimal focusing inside a random sample, (a) Intensity distribution in the trans¬ 
verse dimension for optimal focusing at (LI 4,0) and (LI 2,0). (b) Spatial variation of the 
contrast fl(x) and eigenchannel participation number M(x). (From Ref. 11081) 


The lower value of contrast for focusing at large depth into the scattering sample is related to 
the increasing mesoscopic intensity correlation towards the output surface of the sample. The 





















focusing in the sample is given by Eq. 5 with K replaced k(x), the degree of long range intensity 
correlation at depth x. We expect that the contrast of focusing inside a sample is given in terms 
of the participation number of eigenvalues of the energy density at a depth x into the sample 
M(x) of the field matrix connecting the field at the incident plane and the focal plane at depth 
x. This hypothesis is confirmed in the plot of spatial variation of the focusing contrast fj. (x) and 
1/(1/M(x) — 1 /N) in Fig. 13(b). Good agreement is seen once the depth into the sample is 
beyond the transport mean free path so that the field is completely randomized. 

When a short pulse is incident on a scattering medium, the pulse is stretched in time as 
it is scrambled in space. Focusing of a pulse through a random medium was first achieved 


via the time reversal in acoustics [ 1131. Recently, focusing of a pulse in transmission prop¬ 


agating through an opaque sample has been achieved by adjusting the incident wavefront 
to enhance transmission at a focal spot at a selected time delay ]143|jl46] . The tempo¬ 
ral response of an open random medium to an incident pulse can be naturally described in 
terms of the quasi-normal-modes of the system 169 147} 148| . These modes are character¬ 
ized by their central frequencies (O m and linewidths T m , and a Forentzian shape in frequency, 
r m /2/(T m /2 + i(co — ©,„)). Once excited, each of the modes will release its energy at a rate 
determined by its linewidth T m , exp(— T m t). Spectra of the transmission coefficient between 
incident channel a and outgoing channel b can be expressed as a summation of the modes, 


thn(CO) = Y tZ 


r m /2 


( 6 ) 


Here, the matrix ?!” describes the spatial coupling into the m ,h mode via incident point a 
and coupling out to the output point b from the m ,h mode. We emphasize that, in con¬ 
trast to the eigenchannels, the elements tZ of two different modes can be strongly corre¬ 
lated J69) . The Fourier transform of the spectra of the TM yields a time-dependent TM, 
which characterizes the spatial-temporal response of the random medium to an incoming pulse, 
t(t') = 'Lm t 'b a ex P( — iObJ') ex p{~ r m f'/2), where t' indicates the time delay. To have a pulse 
emerge in transmission at an output point /3 at a certain time delay t', one simply has to ap¬ 
ply phase conjugation of the TM at time delay t', tp a (t'), so that different pulses from incident 
channels a will be in phase at target point fi at time delay t'. 

The spatial-temporal control of pulse transmission through a scattering medium via phase 
conjugation of the time-dependent TM is presented in Fig. 14. The intensity that would be 
delivered to a point at the center of the output surface of the waveguide /3=(0,0) at a selected 
time t' if the transmission matrix were phase conjugated at time t' = 33 ns is presented in Figs. 
14b. A sharp pulse emerges at the selected time delay with intensity peaked at /3=(0,0). The 
temporal profile of the focused pulse is the square of the field correlation function in time, 
|T/(Af)| 2 , where Ff? = (E a (t')E^(t' + + A?))) 1 / 2 . For an incident Gaussian 

pulse, the square of the field correlation function is equal to the intensity profile of the incident 
pulse and is independent of delay time. 

The spatial contrast of the focusing of the pulse is given by 

= (7) 

where M' is the eigenchannel participation number of the TM with a size of N' at time delay t'. 
The time evolution of (M r ) and (j-i) are shown in Fig. 15. Near the arrival time of the ballistic 
wave, the value of M' is close to unity and the contrast is not described by Eq. 7. This is 
because ballistic wave is associated with the propagating waveguides modes with the highest 
speed and therefore the transmitted field is not randomized. Once the transmitted wave at the 
output is multiply scattered, a random speckle pattern develops and the measured contrast is in 
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Fig. 14. Controlling pulse transmission through a random medium with a time dependent 
TM. Spectra of the TM for microwave radiation propagating through a random waveguide 
with a length L = 61 cm are measured from 14.7 to 15.7 GHz in 3200 steps. The TM is 
determined with a single polarization between pairs of 45 points on the incident and output 
surfaces. The time dependent transmission matrix is obtained from spectra of the transmit¬ 
ted field between all points a and b, tb a (v). These spectra are multiplied by a Gaussian pulse 
centered in the measured spectrum at Vo =15.2 GHz with bandwidth c v =150 MHz and 
then Fourier transformed into the time domain. This gives the time response at the detector 
to an incident Gaussian pulse launched by a source antenna with bandwidth a t =\/2n/a v . 
(a) The time variation I bait 1 ) of an incident pulse launched at the center of the input surface 
and detected at the center of the output surface in a single realization of the random sample 
and the average of the time-of-flight distribution {/(?')}. O’) Phase conjugation is applied 
numerically to the same configuration as in (a) to focus at t'= 33 ns at the center of the 
output surface. The Whittaker-Shannon sampling theorem is used to obtain high-resolution 
spatial intensity patterns shown in the inset of (bj. (From Ref. |146|) 


accord with Eq. 7. After the arrival of ballistic waves, the value of M'{t') is seen in Fig. 15 to 
increase rapidly before falling slowly. This reflects the distribution of lifetimes and the degree of 
correlation in the speckle patterns of quasi-normal modes J67]|68) . Just after the ballistic pulse, 
transmission is dominated by the shortest-lived modes [148-150). These modes are especially 
short lived and strongly transmitting because they are extended across the sample as a result of 
coupling between resonant centers. Sets of extended modes that are close in frequency could be 
expected to have similar speckle patterns in transmission, so that a number of such modes might 
then contribute to a single transmission channel (69) , As a result, the number of independent 
eigenchannels of the transmission matrix contributing substantially to transmission would be 
relatively small at early times and M' would be low. At late times, only the long-lived modes 
contribute appreciably to the transmission [i 1511 and so the value of M' will be close to 1. 
Thus for intermediate times, modes with wider distribution of lifetimes than at either early or 
late times contribute to transmission and these modes are less strongly correlated than at early 
times so that M' and the contrast are peaked. 


4. Dwell time and densities of states of the transmission eigenchannels 

Spectra measurements of the TM allow us to explore the photon dwell time of the transmis¬ 
sion eigenchannels. The delay time of transmission in ID disordered samples is given by 
f I{t)tdt / f I(t)dt 1152]. In the limit of vanishing bandwidth, it is equal to the derivative of 
the phase accumulated as the wave propagates through the sample with respect to the angular 
frequency ft), (f)' = d&/d(0 \ 152 1 . Because the delay time in transmission equals the delay time 
in reflection, (j)' is also the DOS divided by n at angular frequency ft), p (ft)). For a multi-channel 
system, (f>! is still the delay time between input channel a and output channel b, but it is not 
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Fig. 15. Time evolution of the maximal focusing contrast (/i). g is well described by Eq. 7 
after the time of the ballistic arrival, t' ~ 21 ns. At early times, the signal-to-noise ratio is 
too low to analyze the transmission matrix. (From Ref. p46| ) 


loner the DOS. This can be seen by noting that a phase shift of order unity occurs within the 
field correlation frequency, in which typically a single peak occurs in transmission, but approx¬ 
imately 8 modes contribute to this single peak. Thus the number of modes is much greater than 
the phase shift. Krein, Birman, Yafaev and Schwinger 1153 155) have shown that the DOS can 
be expressed in terms of the scattering matrix S of the system, p(ft>) = —i / (2k)T r{S^dS / 'dco), 
in which —iS^dS/dco is known as the Wigner-Smith delay-time matrix [156 1571. Measure¬ 
ments of the DOS should therefore involve measurements of the complete scattering matrix. 
For a system with time-reversal symmetry, Brandbyge and Tsukada 1 158) demonstrated that the 
local DOS in an electronic system can be found from the composite derivative of the phase of 
the transmitted field with shift in the potential at a point inside the sample. We found that [jl 15[, 


P(®) = 


£de„ 

~ dco ' 


( 8 ) 


Here, d9„/dco = — i(u^du n /dco — \^d\ n /dco) is eigenchannel density of states (EDOS), which 
is the contribution to the DOS of the n th eigenchannel and is the dwell time of the n th eigen¬ 
channel. d9 n /dco is thus a complementary set of parameters to T„. Thus spectra of the TM yield 
the dynamics of transmission as well as static transmission for each eigenchannel. 

Spectra of z n and d9 n /dco for a diffusive sample with g = 6.9 are presented in Fig. 16. The 
average value of d9 n /dco is seen to decrease as the eigenchannel channel index n increases, 
indicating that the average dwell time becomes shorter as the value of transmission decreases. 
The sum of d9 n /dco represents the DOS of the sample, when the impact of incompleteness of 
the TM is low. We therefore expect that Eq. 8 gives the DOS in samples for which N' M, as 
is readily achieved in samples in which the wave is localized. The DOS of the sample can also 
be found by identifying and counting the contribution of the each quasi-normal mode inside the 
sample. For an open system, the DOS can be expressed as, 


p(®) = La»(®) 

m 


i r „,/2 

K (r m /2) 2 + (to - CO m ) 2 ' 


(9) 


The central frequencies and linewidths of the modes have been found for localized samples 
where the modal overlap is small, by decomposing the transmitted field speckle pattern into a 
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Fig. 16. Spectra of transmission eigenvalues z n and the eigenchannel dwell times dO n /d(Q 
for n = 1,5, 15, 25 for a diffusive sample with g = 6.9. The dwell time for the eigenchannels 
is the contribution of the eigenchannel to the DOS. (From Ref. m3) 


sum of the same set of modes [[69]. This allows us to directly measure the DOS of the disordered 
sample. Good agreement based upon these two approaches for the localized sample with g = 
0.37 is seen and demonstrated in Fig. 17. 




Fig. 17. Comparison of two approaches for finding the DOS of a disordered material, (a) 
Contributions of the individual modes in Eq. 9 to the DOS. The integral of each mode over 
the angular frequency is unity, (b) DOS determined from the TM (red curve) by summing 
spectra of dd n /dco and modes in (a). (From Ref. [115)) 


Since the photon dwell time and the EDOS of an eigenchannel are proportional to the energy 
stored inside the sample 11591, this suggests that the energy stored inside is smaller for low 
transmitting eigenchannels. This is consistent with the results found simulations by Choi el. 
al., EE) and measurements by Gerardin el. al. [91] in single samples. An expression of the 
profile of an eigenchannel within a random medium would be of importance, since it would 
afford a universal description for wave propagation that encompasses the wave inside the sam¬ 
ple as well as the transmitted wave. This would also provide preliminary sense of the intensity 
inside individual eigenchannels channels for a host of application involving energy deposition 
inside random media. This may also be exploited to lower the lasing threshold of a diffusive 
random laser [ jl60||l6l) , in which the threshold is particularly high because the pump energy 
cannot penetrate deep into the sample ]162| . The profiles of eigenchannels over a wide range of 
transmission obtained in 2D simulations using the recursive Green’s function method are shown 


































in Fig. 18. The profiles are averaged over a collection of samples with the same transmission 
T. The profile of the complete transmission eigenchannel with T = 1 is found to be closely re¬ 
lated to the probability of a wave returning to a cross-section at a depth x in an open disordered 
medium 170 116f 118| . We have found that the profile of an eigenchannel with T < 1 can be 
expressed as a product of the profile of the fully transmitting eigenchannel and a function gov¬ 
erned only by the auxiliary localization length which was previously used to parameterize 
the corresponding transmission eigenvalues via T = 1/cosh 2 (L/^'). These results not only give 
a physical meaning to the set of but also enable control of the energy deposition inside the 
opaque media. 



Fig. 18. Spatial profile of the energy distribution of transmission eigenchannel inside the 
sample, (a) Ensemble averages of the eigenchannel energy density profiles W x (x) for eigen¬ 
values 1 = 1, 0.5,0.1, and 0.001 for a diffusive sample withL/<^ = 0.05. W x {x) is the energy 
integrated over the transverse dimension and is normalized to equal T at the output surface, 
(b) Ensemble averages of energy density profiles for all transmission eigenchannels with 
the eigenchannel indices n from 1 to N. The linear falloff of the average of the energy 
density over all eigenchannels is in accord with the diffusion theory. (From Ref. |118|) 


5. Conclusions 

In this review, we have discussed the power of the transmission matrix to describe the statistics 
of transmission and control of waves propagating through disordered systems. We have shown 
that the joint probability distribution P(T,M ) determines the statistics of transmission over a 
random ensemble up to second order, while the statistics within a single sample of given T de¬ 
pends only on M. The statistics of transmission in multi-channel disordered systems approach 
predictions for ID, when the transmission is dominated by a single transmission eigenchan¬ 
nel and so M —> 1. In this limit, the SPS hypothesis for scaling in random ID system holds 
for multi-channel disordered samples. The contrast of optimal focusing through and within a 
scattering sample is equal to M of the measured TM, provided that the number of the meas¬ 
ured channels of the TM is much greater than g. The spatial profile of the focused wave is 
given in terms of the degree of long-range intensity correlation and the spatial field correlation 
function while the temporal profile of a focused pulse is equal to the square of the field-field 
correlation function in time. In addition to the transmission eigenvalues, which describe static 
transmission, there exists a complementary set of parameters, which gives the dwell time and 
the integrated intensity inside the sample for transmission eigenchannels, and the contribution 
of each eigenchannel to the DOS. The TM approach allows us to determine the DOS of a 
disordered material without finding each of the QNMs in the system, which can be difficult 
for diffusive samples because of the modal overlap. The ability to deposit energy deep within 










the sample via high-transmission eigenchannels may find applications in deep tissue imaging, 
depth profiling, photodynamic therapy and random lasing. 
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